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Abstract — Two-level system fluctuators (TLS's) in the tunnel 
barrier of a Josephson junction have recently been demonstrated 
to cause novel energy splittings in spectroscopic measurements 
of superconducting phase qubits. With their strong coupling to 
the Josephson junction and relatively long decoherence times, 
TLS's can be considered as potential qubits and demonstrate 
coherent quantum effects. Here, we study the effective interaction 
between the TLS qubits that is mediated by a Josephson junction 
resonator driven by an external microwave source. This effective 
interaction can enable controlled quantum logic gates between 
the TLS's. Our study can be extended to other superconducting 
resonators coupling with TLS's. 

Index Terms — Superconducting resonators, Quantum theory, 
Superconducting device noise 



I. Introduction 

Two-level system fluctuators (TLS's) are a ubiquitous source 
of decoherence for solid-state qubits 0, 0. In supercon- 
ducting qubits 0, TLS's have been widely studied both 
experimentally and theoretically and are often considered as 
the source of low-frequency (1//) noise 0, 0, 0, 0, 0, 
0, IflOll . In recent experiments, energy splittings in spectro- 
scopic measurements have been observed in superconducting 
qubits, showing coherent coupling between the TLS's and the 
qubits ifTTl . fi2l . The TLS's have demonstrated decoherence 
times much longer than that of the superconducting qubits 1 13 1 
and hence can themselves be considered as effective qubits for 
testing quantum information protocols IfTlll . |[T21l . fBl . |fT31l . 

ma, ini, ma, ma, ma. 

Josephson junctions can be operated as microwave res- 
onators no, (si, m, mg, eh, 00, 02, ma, 1291, ma, 

(3T). In a previous work 11321 , a cavity QED approach ||33l to 
characterizing the coupling between TLS's inside a Josephson 
junction and the junction resonator was suggested where a 
Jaynes-Cummings model was derived and the coupling be- 
tween the TLS's and the resonator can be modulated with an 
applied magnetic field. By measuring microwave transmissions 
in the junction resonator, various properties of the TLS's can 
be probed including their spatial distribution and the coupling 
mechanisms to the junction. 

In the following, we study the effective interaction between 
the TLS's that is mediated by the junction resonator. In our 
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system, because the junction resonator can have a decay rate 
that is much stronger than that of the TLS's, the effect of 
the resonator decay needs to be taken into account. We will 
present the effective interaction both for coupling with a 
high-Q resonator and for coupling with a strongly damped- 
resonator. A microwave source can be applied to the junction 
resonator [32] which provides us with a tool to control the 
coupling between the TLS's and the resonator. This paper is 
organized as follows. In Sec. II, we present the model and of 
the coupling between the TLS's and the junction resonator . In 
Sec. Ill, we derive the effective interaction between the TLS's 
mediated by the resonator mode. In Sec. IV, we will briefly 
discuss the implementation of quantum logic gates between 
the TLS qubits and the decoherence effect. We discuss the 
readout of the TLS qubits in Sec. V and give the conclusions 
in Sec. VI. 

II. Model 

The system is depicted in Fig. Q] A Josephson junction can 
be described in terms of the gauge invariance phase <£> and 
its conjugate momentum P$ ||34l with a capacitive energy 
Pj/2Co and a potential energy — Ej cos(2e$/ft), where Co 
is the total capacitance and Ej is the Josephson energy. When 
combined with the inductance L in an RF SQUID loop, the 
Hamiltonian of the resonator can be written as 



R c = ^k- E J cos(2e$/ft) 



2L 



(1) 



with external magnetic flux $ ex inside the SQUID loop. The 
Hamiltonian describes an oscillator mode: H c w P|/(2Co) + 
C uj 2 ($ — <J> s ) 2 /2 with a shift $ s and a frequency 



4e 2 Ej cos(2e$ s /n) 



LC 



h 2 C 



(2) 



both depending on the magnetic flux $ el . 

Different coupling mechanisms between TLS's and the junc- 
tion have been discussed and observed. For example, TLS's 
can couple with the critical current of the junction in the form 
— {2e/h)Ej$J2„ jn ' where j n and a n are, respectively, 
the polarization vector and the vector of the Pauli spin matrices 
for the n th TLS. Below we assume j n — (j x ,0,0) for 
simplicity. To the lowest order of $ — $ s and with the notation 
$ - $ s = ^h/(2C Lu c )(a + at), the total Hamiltonian of the 
coupled system is 



Hi 



fvaj c a) a + Hi + e(t)(a + a') 



Ehio n 



+ gn(acr n + + a T cr„_) 



(3) 
(4) 
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Fig. 1. Two-level systems (TLS's) inside a Josephson junction. The dashed 
arrows are the magnetic field. The letter 'i' indicates the insulating layer and 
V the superconducting layers. 



where a (aJ) is the annihilation (creation) operator of the 
resonator mode, u n is the frequency of the rfi 1 TLS, and 
<7n = Ejj xn ^Jh/ \2Cquj c ) sin(2eQ s /h) is the coupling con- 
stant. The Hamiltonian includes a driving on the resonator 
mode with a driving amplitude e(t). Eq. (O has the form of a 
Jaynes-Cummings model that has been widely studied in cav- 
ity QED systems. Note that coupling with electrical (dielectric) 
field inside the junction can be derived similarly ifTTl . [32]. 

Dissipative effects can be an important factor when we study 
the effective coupling between TLS's. The decay rate of a 
junction resonator is much stronger than that of the TLS's. 
In the following, we treat the resonator decay as a bosonic 
bath [35] that couples to the resonator with the Hamiltonian 
H K = Yl f^>ka\a k + c k {a\a + a^a k ), where a k (a\) is the 
annihilation (creation) operator for the bath modes, uo k is the 
frequency of the modes, and c k is the coupling constant that 
is related to the decay rate k with: 7r^^ S(u> — uj k ) = k. In 
the master equation approach, the decay can be described in 
the Lindblad forms as 



Op 

— = -i[H t , p] + K,£(a)p. 



(5) 



In this work, we are interested in TLS's that are in the same 
frequency range as that of the resonator, i.e. gegahertz, so 
thermal fluctuations can be neglected. The Lindblad forms can 
hence be expressed as C(o) = i(2o£ot — Co^ o— oC) for an 
operator o. The intrinsic noise bath of the TLS's is neglected 
given their long decoherence times. 

III. Effective Interaction 

TLS's do not interact directly due to their low density inside 
the Josephson junction. However, because of their coupling 
to the same cavity mode, an effective interaction can be 
obtained. In this paper, we study the effective interaction in 
two situations: TLS's coupling with a high-Q resonator with 
moderate to high decay rate ll36l . and TLS's coupling with 
a strongly-damped resonator where the resonator decay is 
stronger than the detuning and the coupling between the TLS's 
and the resonator. 

High-Q resonator. We consider the dispersive regime |37|, 
1 3 8 1 , |39| where the resonator is far detuned from the TLS's 
with g n <C |A nc |. Here, A„ c = A„ — A c is the detuning 
between the n th TLS and the resonator mode. The effective 
interaction can be derived by applying a unitary transformation 



to the Hamiltonian: H t = UH t U^ where the transformation 
is 

U — e - f ( a - at )/ A c J~J e -9„(aV„_-CT„ + a)/A no _ ^ 



The total Hamiltonian then becomes H t = H c + H e ff + H x , 
where H c = hA c a^a denotes the Hamiltonian of the resonator, 
H e ff the effective Hamiltonian of the TLS's, and H x the 
residual coupling between the TLS's and the resonator ll36ll . 
We derive 



Heff=Yl 



hA n 



a, 



H k (7) 



which includes single qubit terms with effective detuning 
A„ = A„ + (gl/A nc )(l - 2e/A c ) and effective Rabi fre- 
quency fl nx — 2eg n /A nc , an effective interaction Hi nt = 
^mni&n+Cm- + cr m+ cr„_)/2 with the coupling constant 



A, 



2A A ' 



(8) 



and an induced coupling to the bath modes of the resonator 
Hk = Yjn,k(9nCkl A nc ){a n+ a k +a\a n -). The residual cou- 
pling has the form 



H r 



E 



9l 



a) a + ' 



A c - 2A r 



(9) 



with a Stark shift for the resonator and an extra coupling to 
the resonator amplitude. The unitary transformation shifts the 
amplitude of the resonator to zero ((a) » 0) for finite driving 
e, so that the effect of the Stark shift is always small. 

Strongly-damped resonator. The decay rate can reach ten's 
of megahertz in lossy resonators and cannot be neglected when 
compared with the coupling between TLS's and the resonator. 
We set the driving to be e(t) = 2e cos ui^t with a frequency ujd 
and an amplitude erj. To study the dynamics in the Heisenberg 
picture, we start from Eq. (0. With 6 = i[H t + H K , o] for an 
arbitrary operator o, we have 

d = -iuj c a - i 2J g n o~n- - ie - i c k a k (10) 
a k = -ic k a - iuj k a k (11) 

which gives 

a k = a k (0)e-^ kt - ic k J dt 'e"^-*') a(t'), (12) 

with a k (0) being the noise operator in the Schodinger picture. 
Substituting Eq. ( fT2l into Eq. ( fTOb . we derive the following 
relation for the resonator mode: 



a = — iA r 



ie 



it/kg 



(13) 



where A r 



lo c i is the detuning of the resonator and 



(lMF)£o*(0)« 



-iujkt 



is the input field of the bath 



modes. Similar equations can be derived for the TLS's: 

<j n _ = -iA„o- n _ + ig n a mz a 
<? n + = +iA n a n+ - ig n a'a mz 
Onz = 2ig n a'o- n - - 2ig n a n+ a 



(14) 
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where A„ = uj n —LOd is the detuning of the n th TLS. In the bad 
cavity limit with n ~ A nc ,g n 3> J n iiln2, we can eliminate 
the resonator mode by setting the right hand side of Eq. (|T3T > 
to zero. This gives 



(15) 



where the resonator adiabatically follows the dynamics of the 
TLS's. The conjugate relation for a> can be derived as well. 
Now substituting Eq. (15[ into Eq. (fl4] i. we can derive a set 
of equations that govern the dynamics of the TLS's: 




B„(16) 



where A n determines the dynamics of a single TLS, i.e. 
determines the parameters in the Bloch Equation for a single 
TLS, with 



-iA n - 72 iSl n + A„ 

A„ = | iA* n - 72 -ifi* + At 

2zQ* - 2At -2ifi n - 2A„ - 7l 

and £>„ determines the effective interaction with 

2zA mrl (J m _|_f7 n _ 2zA nm C7 n _|_(7 m _ 

The parameters for a single TLS in matrix A n are: the effec- 
tive detuning A„ = A„ — A c g 2 /(k 2 + A 2 ), the effective Rabi 
frequency il n = —ig n eo / (n+iA c ), the induced dephasing rate 
72 = <? 2 k/(k 2 + A 2 .), and an induced decay rate 71 = 272. 
The induced dephasing (decay) is due to the effective bath 



A, 



fjn 



iA r 



(17) 



in matrix A n . For TLS's to exhibit quantum coherence, the 
decoherence rates 712 must be weaker than the other time- 
scales in the system, as will be discussed below. The effective 
coupling constant can be derived from B n with 



A, 



K + *A C 



(18) 



and satisfies \ nm = A^„. The coupling depends on A c in a 
similar way as does the Rabi frequency Vt n . 

The following effective Hamiltonian can then be derived for 
the TLS's: 



H, 



eff - _ 

+ y ^nmVn+Vm- + A^ m (T m +<7„_ 

(n,m) 



(19) 



An interesting difference between the coupling in Eq. ((SJ and 
in Eq. ( fT8l is that the effective coupling in Eq. (|T8b doesn't 
depend on the frequencies of the TLS's. In Figj2] we plot the 
magnitude of these two couplings for comparison. 
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Fig. 2. Couplings versus detuning. Dash-dotted line: coupling in Eq. (8) and 
solid line: coupling in Eq. {18}. The parameters are — 2vr X 10 MHz, 
A m = 2tt X 32 MHz, and k = 2tt x 5.1MHz. 



IV. Quantum Logic Operations 

The above system can be used to implement universal 
quantum operations on the TLS's. The driving amplitude, 
driving frequency, and the resonator frequency oj c (or detuning 
A c ) can be independently adjusted. Choosing two linearly- 
independent Hamiltonians Hi and H2 from the general expres- 
sion of H e ff, a complete set of operators can be constructed 
from the commutators such as [Hi,H2\, [Hi,[Hi, H2]], 
etc. [40|. This shows that universal quantum gates can be 
realized by adjusting the above parameters. Details on how to 
realized the gates such as the SWAP gate and the Hadamard 
gate can be found elsewhere [36|. 

The coupling between the TLS's and the resonator induces 
extra decoherence on the TLS's due to the decay of the 
junction resonator. For coupling with a high-Q resonator, a 
noise term H K is generated; for coupling with a strongly- 
damped resonator, a term A„ is generated. To realize quantum 
logic gates, it is necessary that the decoherence rates are much 
smaller than the effective Rabi frequency £l nx and the effective 
coupling constant A mn . It can be shown that the decoherence 
rate for coupling with a high-Q resonator is t7 ~ g^nj A 2 nc - 
With k |A nc |, hundreds of quantum operations can be 
performed within the decoherence time even when the decay 
rate of a resonator is a few megahertz, as has been studied in 
detail in Ref . Il36l . However, the decoherence rate for coupling 
with a strongly-damped resonator is t^ 1 ~ gf l n/(n 2 + A 2 ). 
With k > A c , we have t^ 1 ~ fl n x,^mn, and only a few 
operations can be performed at the best. 



V. Readout 

The junction resonator can function as a readout device for 
the TLS's. We consider the measurement of the n th TLS in 
the dispersive regime. Let the frequency of the resonator be 
close to the frequency of this TLS, but have an off-resonance 
that satisfies g n |A nc |. All the other TLS's are very far 
off-resonance from the resonator, and so have much smaller 
Stark shifts. This can be achieved by switching the resonator 
frequency and the driving frequency at a nanosecond time- 
scale. Because the dynamics of the TLS's happens at a much 
slower rate than this switching rate, the state of the TLS can 
be considered unaffected during the switching. A measurement 



4 



of the transmission or reflection in the junction resonator can 
be used to reveal the qubit states. 

Meanwhile, phase sensitive detection of the stationary state 
of the resonator can give direct measurement of the TLS's 
for a strongly-damped resonator, according to Eq. ( fT3T >. For 
example, we have 

f _ — 2e A e J2 n 9nKa ny + g n A c a nx 

2 i A2 ^2 i A2 ■ 



a + a' — 



When the couplings (<?„'s) are different for different TLS's, the 
output of the resonator provides a readout of multiple TLS's 
in a single measurement ETI . In addition, by choosing the 
phase of the measured canonical variable of the resonator, we 
can choose which TLS operator will be measured. With <fi = 
a.Tg(—igi/(n + iA c )), we have ae~ 1 ^ + ci^e 1 ^ oc u\ x and a 
measurement of <j\ x is performed. 

This scheme also provides a measurement of the time 
dependence of the properties of the TLS's. It can be shown 
that 

{a){t)a{t)) = TT Str [e iStt a^e- iStt aW(Q)] (21) 

where W(Q) is the initial density matrix including the envi- 
ronmental degrees of freedom and H t is the total Hamiltonian 
including the system and the bath. The trace is taken over 
the system modes (index s) and the bath modes (index r). 
Substituting the expression in Eq. (TT~5T > into Eq. (f2TT > and 
considering one TLS (n — 1) for simplicity, we derive that 

K + A~ 

where C(t) = (ax+(t)ai-(t)) and M(t) = (a 1+ (t)+ai-(t)). 
The measured results for the correlation function C(t) can 
be used to interpret the parameters in the matrix A n to 
reveal the properties of the TLS's via the quantum regression 
theorem [35]. 
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VI. Conclusions 

To conclude, we studied the effective coupling between 
the TLS's inside a Josephson junction using quantum optics 
approaches. Two situations are studied and compared: TLS's 
coupling with a high-Q resonator mode and TLS's coupling 
with a strongly-damped resonator mode. Our results indicate 
that the couplings in these two regime have very different 
properties. Universal quantum gates can be realized on the 
TLS's when they are coupled with a high-Q resonator. While 
the fast decay of a strongly-damped resonator can destroy the 
coherence of the qubits and affect the successful realization of 
the quantum gates. We also discussed the readout of the TLS 
qubits via the detection of the junction resonator. 
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